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A quantitative analysis of the process of condensation of bosons both in harmonic traps and in 
gases is made resorting to two ingredients only: Bose classical distribution and spectral discretness. 
It is shown that in order to take properly into account statistical correlations, temperature must 
be defined from first principles, based on Shannon entropy, and turns out to be equal to f3~^ only 
for T > Tc where the usual results are recovered. Below a new critical temperature Td is found, 
where the specific heat exhibits a sharp spike, similar to the A-peak of superfluidity. 
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In the analysis of low temperature behaviour for Bose 
systems, temperature is typically thought of as given by 
T = l/ksP, (3 being the Lagrange multiplier one has to 
introduce in the customary ensemble approach to account 
for energy conservation. Such definition however im- 
plies in fact the assumption of statistical independence of 
modes, which not only is not necessary but may even be 
in contradiction with the basic features of the conceptual 
scheme leading to the construction of quantum statistics. 
Indeed, the latter rests on the notion of density matrix, 
which is correctly introduced only after the Fock space 
has been properly generated over the tensorized product 
of the appropriate number of single-particle states. This, 
in turn, implies the existence of mode correlations which 
are there even in the absence of interactions. Only when 
temperature is high enough as to break such correlations, 
the Fock density matrix can be lifted to a product of sin- 
gle particle density matrices, in which case, however, the 
only statistics achievable is MaxwcU-Boltzmann's, where 
correlations are absent. 

Ever since 1907 Einstein himself discussing the be- 
havior at low temperature of the quantized harmonic os- 
cillator in the frame of Boltzmann statistics, pointed out 
how classical results on equipartition of energy are, in 
agreement with the Nernst theorem, no longer valid for 
T < hv/ks, by obtaining a good fit of the heat capacity 
data for diamond. One of his basic assumptions was the 
extensivity hypothesis. The latter is usually transferred, 
adopting a tensor product structure for the Hilbert space 
of states of a multi-particle system, to the whole Quan- 
tum Statistical Mechanics (QSM) in spite of the fact that 
it can be justified for the Gibbs-Boltzmann distribution 
only, related in this case to the statistical independence 
of different modes. It worked well for Einstein when he 
dealt with the carbon atoms in diamond (which of course 
do not condense) in that, being fixed at their positions 
in the crystal, these are distinguishable, while one of the 
characteristic features of quantum mechanics is indistin- 
guishability of the particles. 

Several authors have discussed Bose Einstein conden- 
sation (BEG) in magnetic traps, all on the basis of that 



same definition of temperature finding - to- 

gether with Tc oc N^^^ - a heat capacity monotonically 
decreasing to zero for T below Tc, in agreement with what 
happens for true gases under the same hypotheses. 

On the other hand Bose, when he "invented" bosons, 
broke scale invariance in the same way as Democritus did 
in the 5th century BC: because of correlations, a conden- 
sate is not the disjoint union of smaller condensates, in 
the same way as an atom (ancient or modern) is not 
equal to two half-atoms. The properties of a condensate 
do depend on the number N of atoms it contains, as one 
cannot consider the thermodynamic limit § 0, 0] 
that would require, together with N — > oo, hv oc N~^'^ 
(i.e. the continuum), in confiict with the properties of 
magnetic traps where the spectrum (discrete) has a spac- 
ing which is fixed and comparable with ksTc. In other 
words the distinction of physical quantities in intensive 
and extensive is not consistent with Bose-Einstein con- 
densation, as proved, among others, by the mentioned 
relation Tc oc N^/^. 

All these well known features are little relevant for 
T > Tc, because the particles there are distributed over 
all levels, and the continuum approximation appears to 
be sufficient both for describing the regime for T > Tc 
and for determining Tc itself (minor corrections may 
arise, that will be considered elsewhere [^). However, as 
soon as the condensate fraction c = na/N becomes rel- 
evant (T < Tc), statistics becomes significantly different 
from Boltzmann's and the very definition of temperature 
must be reconsidered together with, but independently 
from, the question of equipartition. Temperature can be 
defined, in a quite general and basic way, in terms of two 
fundamental physical quantities; internal energy E and 
entropy S. Thus instead of connecting T, because of ex- 
tensivity, to the Lagrange multiplier /3, one can derive it 
from the separate measurements of E and S. 

We shall then consider ab initio, in the grand-canonical 
ensemble, a system of N bosons of total energy E and 
we shall explicitly evaluate the Lagrange multipliers as 
functions of N and E. Successively, Shannon's form of 
entropy will be our starting point to derive all equilibrium 
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thermodynamic properties. In other words, our hypothe- 
ses will be the three consistent assumptions: 
a) The Bose-Einstein statistics: 



b) a discrete energy spectrum (as suggested by the typical 
experimental set-up with magnetic traps, where ksTc ~ 
10 hi'), and 

c) Shannon entropy, evaluated in terms of N and E. 
Temperature T = T{N, E) will then be obtained as a 
derived quantity by T = dE/dS. 
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Fig. 1. Specific lieat C/fes vs. T/T^ for N = 10® bosons in an 
isotropic harmonic trap, as predicted by the entropy approach. 

The resulting thermodynamics of bosons at tempera- 
tures above and below Tc is thus simply derived numer- 
ically, both for an harmonic potential and a cubic box, 
from the two inputs N and E. The temperature scale Tc 
is found to confirm the results of reference ^ , the specific 
heat C/kB has the appropriate value for T > T^, 3 for 
the three dimensional harmonic oscillator and 3/2 for the 
gas in a box (as it should, since the system does not con- 
dense and behaves essentially as a classical gas) ; however 
new features appear for < T < Tc. For both the phys- 
ical situations considered, heat capacity exhibits an un- 
expected behaviour; at Tc, it falls briskly down to almost 
zero (entropy and energy decrease with quite different 
slopes because only a few atoms migrate progressively 
from the warmer tail of the spectrum into the ground 
level) and raises up abruptly at a new critical tempera- 
ture Td (where uq/N « 1/2), at which it exhibits a sharp 
spike, analogous to that of superfluids As no in- 

teraction has been introduced, other than the confining 
potential, this enhance the idea that superfluid transition 
may indeed be ascribable, as originally suggested by On- 
sager and Penrose essentially to Bose-Einstein con- 
densation. For T < Td, specific heat goes to zero when 
T goes to zero, in agreement with Nernst's theorem. It 
is suggestive that, for T ^ 0, C cx T"^ for harmonic traps 
whereas C oc T^/'^ for bosons in a box. In Fig.l C/ks 



is given vs. the "probabilistic" temperature. To check 
the procedure, C has also been numerically evaluated as 
dE/dp~^, checking that this coincides with the standard 
textbook form (see, e.g., ||l2|). It is thus clear that the 
effect lies in the different definitions of temperature and 
not in computations. 

We shall now summarize a few technicalities of the pro- 
cedure, emphasizing that the scheme proposed to deal 
with correlations could in principle be utilized also for 
interactions. 

1) A set of values for N in the range 10^ 10® as well 
as a set of values for the mean energy per particle E/N 
in the range from the ground state up to 10^ units hi' 
for the magnetically trapped system and /LP' for the 
cubic box (of side L) have been considered. The differ- 
ence between the two cases lies in the energy spectrum 
and its level multiplicity. The full range of condensation 
< c = uq/N < 1 has thus been covered. 

2) For each pair of values N and E the system of two 
equations in the two unknowns a and /3: 

'^rii{a,/3)^N , ^^(q;, (3) = E , 

(in the sums the correct state multiplicities have of course 
been taken into account) was solved. The two functions 
a = a{N, E) and /? = I3{N, E) were thus found. 

3) Successively, the Hi = ni{N,E) were obtained from 
eq. (0). In particular the condensate fraction c = rio/N 
was determined as a function of iV and E. 

4) The probability of finding a particle in state i, pick- 
ing it up at random in the ensemble is Pi{N, E) — rii/N . 
This allows us, adopting Shannon's definition, to find the 
entropy 5 as a function of N and E; 

S = -ks ^ Pi Inpj. 

In this way the correlations induced by quantum statis- 
tics enter into play only in the definition of the 

Pi's which are such that p^"'^'^'* ^ p't^pf''' if the system 
is split into the union of two parts. It is interesting 
to note that, had we considered Boltzmann distribution 
and not eq.(|l]), Gibbs' distribution would have been ob- 
tained and the results of 0| reproduced. Once the piS 
are known, any further analysis of equilibrium thermo- 
dynamical properties is mere statistics. 

5) Temperature T is obtained, as a function of TV and 
E, by numerical derivation: 

r.T,iv.E).[||]r 

6) The relation T ~ T{N, E) is monotonic in E for 
fixed N, and it can be inverted to obtain E — E{N,T), 
whence the specific heat is obtained by (numerical) dif- 
ferentiation (see Fig. 1): 
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c = 



' djE/N) ' 
dT 



N 



7) "Probabilistic" temperature T, despite its concep- 
tual relevance, is perhaps not the physical parameter that 
best describes the process of condensation: it is diffi- 
cult to measure in the experimental set-up of traps and, 
as one does not have equipartition, its intuitive mean- 
ing of "average energy per particle" is lost. On the 
other hand, experimentalists define temperature deriving 
it from the high energy tail of the velocity distribution 
for which they assume equipartition. In other words, 
the "experimental" temperature is 
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3kB N — no 



We therefore describe the onset of condensation by ex- 
hibiting the behavior of c vs. (Fig. 2), that we compare 
with c vs. T (Fig. 3). 
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Fig. 3. c vs. T/Tc (T probabilistic temperature). Apparently sim- 
ilar to Fig. 2, just merging to zero in a smoother way, this is indeed 
deeply different, in that the critical temperature is Td instead of Tc 
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Fig. 2. Condensate fraction c = no/N vs. "experimental" tem- 
perature 0/Tc. The curve is perfectly fitted by c = 1 - (0/Tc) . 

For the sake of completeness, we report also the 
"experimental" specific heat C = d{E/N)/dQ (Fig. 4). 
It should be noticed that since above critical temperature 
equipartition holds also in our case, the two critical tem- 
peratures coincide (8c = ~ .9hv /ksN^') and equal 
that given in ||] . In Fig. 5 the specific heat of the "warm 
tail", = d{3kBQ)/dT is shown vs. T/T^). 

8) As we have the complete description of the sys- 
tem in function of N and E, all physical observables 
can be obtained. For instance the chemical potential, 
H = ^{N, E) = (dE/dN)^, which gives a ^ -(3^. 

In order to make the figures more easily comparable 
with known results in the literature, only those referring 
to a single value of N {N = 10^) and to the harmonic 
confining trap are reported. It should be pointed out 
that, taking into account that cx N^^^, no substantial 
differences appear either for different values of N nor for 
the gas in a box. 



Fig. 4. "Experimental" specific heat C/fcs = d(E/N)/d{kBB) vs. 
S/Tc. Well reproduces the result in 

Our interest is focussed on different aspects. On the 
one side the proposed combinatorial approach empha- 
sizes the relevance, for a proper interpretation of phe- 
nomenological data, of a definition of temperature ade- 
quate to the experimental set up. On the other, it hints 
to the existence at very low temperature of effects that 
would be worth investigating in the laboratory. Such 
effects are probably at present beyond reach, but exper- 
iments are in constant evolution, and the effort might 
open interesting new perspectives. Finally, it may pro- 
vide a novel insight on the A-transition in superfluids. 

QSM has always been at the border between statisti- 
cal mechanics (SM) and probability theory. For the latter 
p3[ SM is related to information theory, where also the 
lexicon is different from that usual one in physics. For 
instance, "Boltzmann statistics " is simply the common 
limit for ^ 1 of the only two existing statistics, Bose 
and Fermi in while is just one possible "physical 
distribution" with unlimited occupation numbers, like 
that of Bose in O, [ITsl. 
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Fig. 5. Specific heat of the warm tail Cw/ks vs. T/Tc. Cw/kB 
is constant (= 3) for T > Tc as expected classically, and it is two 
orders of magnitude smaller for < T < Tc, and ~ 4 times larger 
for < T < Trf. 

In distribution theory all information is contained in 
weights such as those in eq. (P, while the fundamental 
hypothesis of statistical mechanics is statistical indepen- 
dence of the subsystems. In QSM it is the symmetriza- 
tion imposed by particle indistinguishability that one has 
to perform, e.g. moving to Fock space, that breaks such 
independence. Results of physical interest are neverthe- 
less typically in agreement, which is due in most cases 
to the fact that statistical correlations are irrelevant and 
therefore all physical distributions are equivalent. Well 
established exceptions are fermions at low temperature 
and, referring to bosons, the black body, where photons 
are well described by eq.(0). In black body theory, how- 
ever, the situation is very peculiar because photons es- 
sentially cannot interact with each other. 

On the other hand, in view of the continuous improve- 
ments of experimental techniques, extending QSM to de- 
scribe experiments where rii 1, the spectrum is dis- 
crete, the number of particles and modes is large but not 
so large as to allow a continuum description, is certainly 
a challenge that will have to be faced soon. Indeed all 
of this is already happening just for one physical set-up: 
Bose condensation in magnetic traps. 

In the condensate most particles are in the fundamen- 
tal level and, in typical experimental situations, ksTc 
is only of the order of ten times the spectrum spacing 
(Tc « lO^nK, hv/ks ~ IQnK). The usual approxima- 
tion in terms of continuum density with ng obtained by 
difference, appears to be improvable. Moreover, the sta- 
tistical independence of levels - acceptable in the black 
body ~ is questionable when the container is a magnetic 
field and has no kind of thermal exchange with the atoms. 
Thermodynamic equilibrium of each mode is reached be- 
cause not of the interaction of the mode itself with the 
bath, but of the jumping of atoms from one level to an- 
other, as clearly shown by the fast restoring of equilib- 
rium after each cooling. In the case of gas, temperature is 
controlled by the walls and the physics of thermalization 



is different but, if the statistics is that of Bose, statistical 
independence must be excluded. 

Keeping the spectral structure into account provides 
the correct conceptual scheme as opposed to a theory 
where continuum and equipartition are assumed a pri- 
ori. On the other hand, many body effects so far not 
considered should be analysed; however we expect that 
including in the picture interactions in the ground state, 
as described in the Gross-Pitaevsky approach [|l6|, [||, 
may modify at most the energy spectrum, and the ensu- 
ing differences in numerical values should not be able to 
substantially mask the effects discussed. 

It should be pointed out that the description presented, 
based on the assumption of spectral discreteness, leads 
not only to qualitative behaviours, but it determines as 
well physical observables: specific heat, filling order pa- 
rameter, chemical potential, etc., all of which can be ex- 
perimentally measured, related to the characteristics of 
the experiment [N , E, v). To this effect, it should be 
observed that in real experiments the harmonic poten- 
tial is usually anisotropic: extension of our results to 
this situation is straightforward and leads to quite simi- 
lar predictions. 
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